Spin-drag relaxation time in one-dimensional spin-polarized Fermi gases 
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Spin propagation in systems of one-dimensional interacting fermions at finite temperature is in- 
trinsically diffusive. The spreading rate of a spin packet is controlled by a transport coefficient 
termed "spin drag" relaxation time r^d- In this paper we present both numerical and analytical 
calculations of Tsd for a two-component spin-polarized cold Fermi gas trapped inside a tight atomic 
waveguide. At low temperatures we find an activation law for Tad, in agreement with earlier cal- 
culations of Coulomb drag between slightly asymmetric quantum wires, but with a different and 
much stronger temperature dependence of the prefactor. Our results provide a fundamental input 
for microscopic time-dependent spin-density functional theory calculations of spin transport in ID 
inhomogeneous systems of interacting fermions. 
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I. INTRODUCTION 

Quantum many-body systems of one-dimensional (ID) 
interacting particles have attracted an enormous interest 
for more than fifty yearai and are nowadays available 
in a large number of different laboratory systems rang- 
ing from single-wall carbon nanotubes^ to semiconduc- 
tor nanowires^, conducting molecules^, chiral Luttinger 
liquids at fractional quantum Hall edges^, and trapped 
atomic gascsS.. 

Regardless of statistics, the effective low-energy de- 
scription of these systems is based on a harmonic the- 
ory of long- wavelength fluctuations^, i.e. on the "Lut- 
tinger liquid" modeli. The distinctive feature of the Lut- 
tinger liquid is that its low-energy excitations are inde- 
pendent collective oscillations of the charge density or 
the spin density, as opposed to individual quasiparticles 
that carry both charge and spin. This leads immediately 
to the phenomenon of spin-charge separatior^i, i.e. the 
fact that the low-energy spin and charge excitations of 
ID interacting fermions are completely decoupled and 
propagate with different velocities. This behavior has 
been recently demonstrated by KoUath et al^ in a nu- 
merical time-dependent density-matrix renormalization- 
group study of the \D Hubbard model for fermions, by 
Kleine et al^ in a similar study of the two-component 
Bose-Hubbard model, and, analytically, by Kecke et alJ^ 
for interacting fermions in a ID harmonic trap. The pos- 
sibility of studying these phenomena experimentally in 
ID two-component cold Fermi gases^, where "spin" and 
"charge" refer, respectively, to two internal (hyperfine) 
atomic states and to the atomic mass density, was first 
highlighted by Rccati et alr^. 

In a recent paper— two of us have pointed out a 
new aspect of spin-charge separation: namely, spin ex- 
citations are intrinsically damped at finite temperature, 
while charge excitations are not. The physical reason for 
this difference is easy to grasp. In a pure spin pulse the 
up-spin and down-spin components of the current are al- 
ways equal and oppositely directed, so that the density 



remains constant. The relative motion of the two compo- 
nents gives rise to a form of friction known, in electronic 
systems, as "spin Coulomb drag"^^^^^Al^^. This is re- 
sponsible for a hydrodynamic behavior due to a random- 
ization of "spin momentum" deriving from excitations in 
both the g — > and q — *■ 2fcF regions. In the course of 
time this leads to diffusive spreading of the packet. No 
such effect is present in the propagation of pure density 
pulses, which are therefore essentially free of diffusion in 
the long wavelength limit^ (by contrast, a density pulse 
in a normal Fermi liquid is always expected to decay into 
particle-hole pairs - a process known as Landau damp- 
ingi^). 

The analysis of Ref. [H was limited to spin- 
compensated (unpolarized) systems. In particular, the 
calculation of the spin-drag relaxation time Tsd was done 
only for unpolarized systems. It is of great interest to 
extend the calculation to spin-polarized systems. First, 
a finite degree of spin-polarization is the most general 
situation in experiments. Second, and most important, 
any attempt to study the dynamics of spin packets be- 
yond the linear approximation will have to take into 
account the spin-polarization of the liquid within the 
packet. This is particularly crucial in the time-dependent 
spin-current-density functional approachi^, which treats 
spin-current relaxation through a local spin-drag relax- 
ation rate. Calculations done by this method will there- 
fore require knowledge of Tgd function of the local 
spin polarization and temperature. Finally, our study 
is relevant to the closely related problem of the regular 
Coulomb drag (no spin involved) between two semicon- 
ductor wires with different carrier densities. Here the 
two wires play formally the role of the two spin orien- 
tations, and the difference in density corresponds to the 
spin polarization. Our careful analysis of the low tem- 
perature regime yields a new formula for the temperature 
and polarization dependence of Tgd, which is completely 
different from the one that was obtained previously by a 
qualitative method^S. 

The contents of the paper are briefly described as fol- 
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lows. In Sect. M 

we introduce the model Hamiltonian 
and we present our numerical calculations of the spin- 
drag relaxation rate of ID spin-polarized Fermi gases. 
In Sect, mil we report and discuss our main analytical 
results. Finally, in Sect. HVl we summarize our main con- 
clusions. 



II. SPIN-DRAG RELAXATION TIME FOR 
SPIN-POLARIZED SYSTEMS 

We consider a two-component Fermi gas with N atoms 
confined inside a tight atomic waveguide of length L 
(an "atomic quantum wire") along the x direction, re- 
alized e.g. using two overlappingstanding waves along 
the y and the z axis as in Rcf. [g. The atomic waveg- 
uide provides a tight harmonic confinement in the y — z 
plane^i characterized by a large trapping frequency — 
27r X 10 kHz. The two species of fermionic atoms are as- 
sumed to have the same mass m and different spin cr, 
a =1" or J,. The number of fermions with spin ] {N-^) is 
taken to be larger than the number of fermions with spin 
I (iVj), N-\^ > N[. The fermions have quadratic disper- 
sion, Ek = h'^k'^/(2m), and interact via a zero-range s- 
wave potential v{x) = giDS{x). The Fourier transform of 
such a real-space potential is a simple constant, Vq = gn). 
The system is thereby governed by the Yang Hamilto- 
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The effective II? coupling constant 51D can be tuned^ 
by using a magnetic field-induced Feshbach resonance 
to change the 3D scatteri ng length 30. In the limit 
flsD ^ o,±, where a± = \/h?/ (mcu^), one finds giu = 
2fi^ a^D / {ma\)i^. In the thermodynamic limit {Ncr,L 
00, N^/L = Uij) the properties of the system described 
by Ti are determined by the linear density n = n| -|- n^, 
by the degree of spin polarization C, = {N-^ — Ni)/N, 
and by the effective coupling gn). For future purposes 
it will be useful to introduce the dimensionless Yang in- 
teraction parameter 7 = mgn:,/{h^n). We also introduce 
the Fermi wave number k-p = 7rn/2 of the unpolarized 
system, the Fermi velocity vf = hk-p/m, and the Fermi 
energy ep = h'^kp/{2m). 

Within second-order perturbation theory the spin-drag 
relaxation rate (at zero frequency) is given by the for- 
mula— 



~ = / V- '^^9 

Tsd ninimkBl Jq 27r 

^ Jo sinh^[fia;/(2/cBT)] ' 

where StoXo-'''' (9, w) is the spin- resolved finite- 
temperature expression for the imaginary part of 
the ID Lindhard functiorti^, 
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Here kp^ = kp[l + sgn{a)(] [by definition sgn(a-) = +1 
for cr =1 and to —1 for a =|], Na{0) = m/ {nfi^kpa) 
is the spin-resolved density of states at the Fermi level, 
eF[l + sgn((T)C]^, and 



epcr 



mcu q 
v±,a = — — ± 



qhkpcr 2kpa 



(4) 



In Eq. ^ l^a{T) is the spin-resolved chemical potential, 
which is determined by the normalization condition 



dk 



1 



27r exp{[efe - ^i„{T)]/ {kBT)} + 1 



(5) 



In Fig. [T] we report some illustrative numerical re- 
sults for the spin-drag relaxation rate of a ID spin- 
polarized Fermi gas as calculated from Eqs. ©-([S])- Note 
from the inset in Fig. [1] that a finite value of C, has a 
dramatic effect on the low-tempcrature behavior of the 
spin-drag relaxation rate, changing it from linear— to ex- 
ponentially activated. At high temperatures r".^ is seen 



to be insensitive to the degree of spin polarization, be- 
coming asymptotically equal to the unpolarized ^ = 
result. Both observations will be demonstrated in the 
next Section with analytical calculations. 

Note finally that is largest for the unpolarized case, 

simply because the "overlap" between '^mx^^^ {q, to) and 
Qto^x"^ is maximum at = 0. 



III. ANALYTICAL RESULTS 



Low- and high-temperature analytical expressions for 
Tgj^ have been derived in the unpolarized case in Ref. [l^ . 
Here we proceed to derive some analytical results for 
for ^ ^ 0. Similar calculations have been performed 
earlier by Pustilnik et al^ with the aim of studying 
the Coulomb drag between slightly asymmetric quan- 
tum wires. We will comment at length below on the 
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FIG. 1: (Color online) Spin-drag relaxation rate (in units 
of ef/S and divided by 7^) as a function of the reduced tem- 
perature ksT I (2eF) for three different values of the spin polar- 
ization In the inset we show a zoom of the low-temperature 
region < A:Br/(2eF) < 0.5. 



connection between our calculations and the calculations 
reported in Ref. [l^. 

It is crucial to realize^ that a finite degree of spin 
polarization implies the existence of a new tempera- 
ture scale (i.e. in addition to the Fermi temperature 
Tp), that we define as Tc^ = Tp]- — T-pi = 4CTf, where 
Tpo- = £-pa/k^- In what follows we will derive analytical 
results for T ^ nYm{T(^,T-pi). Note that this inequality 
guarantees also that the minority J,-spin component is al- 
ways in the regime of quantum degeneracy, thus allowing 
us to use Eq. ^ also for the imaginary part of the spin- 
resolved density-density response function for the minor- 
ity spin population. 

Before proceeding to illustrate our analytical proce- 
dure [outlined from Eq. ([5]) down below] , it is instructive 
to summarize the main steps followed by Pustilnik et al. 
to derive Eq. (16) in Ref. '2C?. As it will be clear be- 
low, there are two contributions to the inverse spin-drag 
relaxation time, one coming from the low-g region and 
the other one coming from the q « 2kp region, similarly 
to what happens in the unpolarizcd casein. In Ref. |2(J 
the authors focus only on the low-q contribution and use 
the T = rectangular form for the imaginary part of 



the Lindhard response function, i.e. Q'mxa\q, (^)\t=o = 
[■m/{2tfq)]e{hq^/{2m) - |w - vpaql), where Q{x) is the 
Heaviside step function. This immediately implies that 
the product ^■mxf\q,u!)'^mxf\q,LLi) in Eq. ([2]) is it- 
self a rectangle of width Suj = hq^/m — 2i^hqkF /m cen- 
tred at w = qvp. For this reason Pustilnik et al^ 
substitute the thermal factor 1/ sinh^[?ia;/(2fcBT)] with 
1/ siTih^[hqvp/(2kBT)], bringing it outside the uj integra- 
tion in Eq. Then, using the fact that 



a(9,C) 



+ OC 



they arrive at the resulli^S 



duj 'imxf\q,to)^mx''^\q,c^) 
(q - 2CfcF)e(g - 2CfcF) , 



T=0 



(6) 
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Tsd 
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PMGA 2fnT^nkBT j2<:;k 
oc CT exp i-Tc/T) 



dq vf 



q{q - 2CfcF) 

smh^[hqvF/{2kBT)] 
(7) 



for C < 1. 

Using the zero-temperature expression for the product 
'^mxf\q,i^)'^'rnx^l^\q,u>) results, however, in the wrong 

prc-exponcntial factor: the thermal tails of '^mxlr^ {q, <^) 
are exponentially small but they have a finite overlap for 
the two spin populations precisely in the low-frequency 
region where the factor 1/ sinh^ [fko / {2k-BT)] is large and 
thus can originate contributions that are much larger, 
in the low-temperature limit, than the contributions re- 
tained in Eq. Q . The two regions where the thermal tails 

of '^mxf^ {q, i^) and Q-mxf'^ (g, to) overlap are depicted as 
dark grey regions in Fig. [2] and will be discussed at length 
below. In what follows we will show that the functional 
dependence on ( and T of the pre-exponential factor in 
Eq. ([7]) is incorrect. In particular the prefactor T should 
be replaced by the much larger T~^. To show this, we 
proceed as follows. 

Inserting Eq. ^ in Eq. ^ we find, after some straigth- 
forward algebraic manipulations, the following expression 
for ^sd'' 



1 mnl3 f°° , 
— = FTi2 / "-1 ^> 



+ 00 

2 / duj 



I{q,^) 



smh^iPhuj/2) 



Here /3 = (kBT) ^ and I{q, cu) can be written as 



(8) 



w sinh^(/3fi^/2) 

cosli(/36'i)/2 + cosh(/36i2)/2-|-cosli2(/3fia)/2) + 2cosh(/36ii/2)cosh(/36'2/2)cosh(/3fiw/2) 



I 

with We would like to stress at this point that Eq. ^ has 

01 = 4CeF 



muj^ , h^q^ , • (10) 



4 




FIG. 2: (Color online) The one-dimensional particle-hole 
continuum of the majority spin-| component (solid lines) 
and the particle-hole continuum of the minority spin- J, com- 
ponent (dashed lines). The dark grey regions are the in- 
tegration regions that we have considered in our calcula- 
tions. The light grey region is the one considered by Pustil- 
nik et al^. The small thermal-overlap-induced value of 
'Qmx^^\q,i^) X 'Qmx^i\q,ijj) in the dark regions is more than 
compensated by the large value of 1/ smh?[hijj / {2k-BT)]. 



been obtained from Eq. ([2]) without any approxima- 
tion. Two facts are remarkable. To begin with, we 
note that, quite surprisingly, an exact cancellation oc- 
curs between the factor smh^ [pfiw / 2) in the numerator 
of Eq. dH) and the same factor in the denominator of 
the ti;-integrand in Eq. ^ . Secondly, a new temperature 
scale, = 4CeF/fcB, has set in. A color plot of the func- 
tion I{q,uj)/ sinh^(/3?tti;/2), i.e. the integrand in Eq. ([H]), 
is reported in Fig. [3] for two values of C at T = 0.1 Tp. 
We can thus distinguish clearly two regions in the {q, lo) 
plane where this function is non-zero: one region, "close" 
to q ~ 0, which has the shape of an olive leaf and an- 
other region, "close" to 2kp, which has the shape of a 
shark fin. These two regions will be defined accurately 
in what follows. Note also that on decreasing ( the area 
of each of these regions decreases (see also Fig. |3|): in the 
limit C ~* the two regions become two points located at 
(0,0) and (2fcF,0) in the (g, w) plane, as expected from 
the calculations performed in Ref. 

Thanks to the aforementioned cancellation, to make 
some analytical progress in the calculation of Tgd for T — > 
we just need to understand what happens for T — > 
to the four cosh terms in the denominator of Eq. ^ . To 
begin with, note that the first term cosh(/36'i)/2 is inde- 
pendent of q and uj. In the limit T — > we can substitute 
all hyperbolic cosine functions with simple exponentials. 



I{q,uj) T^Q 4 
sinh^(/3ncj/2) exp(/36»i) + exp{P\92\) + exp{(3huj) + exp[/5(6'i 16*21 + nuj)/2] 

4cxp(-/36ii) 

1 + exp[^(|6i2| - 9i)] + exp[l3{huj - Oi)] + exp[/3(|6'2| + hu - 9i)/2] ' 



where we have taken into account the possibility that 62 
becomes negative [contrary to the other arguments of the 
cosh functions in Eq. ^ that arc positive definite] . In the 
limit T — > the exponentials in the denominator of the 
second line of Eq. pT|) are either zero (if their argument is 
negative) or infinity (if their argument is positive). This 
implies that ex-p{/39i)I' {q,iu) differs from zero only when 
the arguments of all the exponentials in the denominator 
of Eq. (fTTj) are negative. Thus I'{q,u!) reduces to a very 
simple form, 

r{q,uj) 4cxp{-i3ei)e{ei-\92\)Q{0i~nuj) 

X 9(^1 - 1^?2| -M- (12) 



More explicitly, I'{q,uj] 
Huj > 0, i.e. 



4exp(-/30i) if 01 - \02\ 



4m 



2£f(1 + C') 



(13) 



and zero elsewhere. Eq. (jl3[) determines two discon- 
nected regions in the (g,w) plane where I'{q,uj) is non- 
zero: (i) a region that we have called above "olive leaf 
(^'ol) enclosed between the two parabolic arches 



^up(g) = 
"^down 



-Sg/h- 

Eq/h- 



- fikpil 



C)q/m 

' C)q/m 



(14) 



with < q < 2C,k-p] and (ii) a region that we have 
called above "shark fin" (Psf) enclosed between the two 
parabolic arches 



^icft('7) = Sq/fi- - hkF{l - C)q/m 



bright {q) 



j/h + hkF{i + C)q/m 



(15) 



with 2(1 — Ofcp < q < 2{1 + Okp. These regions have 
been shown in Fig. [21 vis-a-vis the region implicitly con- 
sidered in Ref. [20. Note that in these domains I'{q,uj) is 
a constant = 4exp(— T^/T), independent of q and uj. 

At this point we can easily calculate analytically the 
asymptotic T — » behavior of the spin-drag relaxation 
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FIG. 3: (Color online) A color density plot of I'{q,ui) as a 
function of q/kp and huj/ep for T = 0.1 Tp. Top panel: ( — 
0.4. Bottom panel: = 0.75. The two regions, olive leaf and 
shark fin, are clearly distinguishable. 
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FIG. 5: (Color online) A comparison between the analyti- 
cal formula in Eq. (|16|) (solid line) and the numerical results 
calculated from Eqs. Q-® (filled squares). 



rate in Eq. ([8]): for Vg = giD wc find 



1 T^o 4m 



g?^ x4cxp(-Tc/T) 
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dqduj 
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exp 
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where we have used that J J, 



I^olUCsf 



dqdoj 



(16) 
(SC^ + 



?>C, /i)k-pe-p /h (sum of the two areas of the olive leaf and 
shark fin domains). In this equation the first contribu- 
tion cx C,"^ comes from the shark fin region (i.e. the re- 
gion "close" to g = 2fcF), while the second contribution 



cx C'^ comes from the olive leaf region (i.e. the region 
"close" to g = 0). The activation law exp(— T^/T) in 
Eq. pB]) is identical to the one reported above in Eq. ([7]) 
(see Ref . [20l) . The pref actor, however, has a completely 
different functional dependence on temperature (propor- 
tional to rather than T) and polarization. Quite 
interestingly, the scaling C ^ T in Eq. pB]) gives the cor- 
rect temperature dependence of the unpolarized result, 
i.e. a linear term from q k, 2k-p and a quadratic term 
from q K Q. In Fig. [5] we show a comparison between 
the analytical formula p6p and the numerical results for 
C = 0.1 as calculated from Eqs. ©-(H]). The agreement 
is clearly excellent in the asymptotic regime T ^ . 

Before concluding, we would like to calculate the high- 
temperature limit of the spin-drag relaxation rate. In this 
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limit the spin-resoved chemical potential is simply given 
b y /Ja.(T oo) = kBT\n[n„XdB{T)], where XdsiT) = 
^J^-kW- I (mkBT) is the thermal de Broglie wavelength. 
Using this result one easily finds 

j_ T^oo m 8 /fcsry^/' 

T,d ^ 2h^n^ 7rV2 2£F / 

X £^d,,.^,e.,^-^), (17) 

which for Vq = giB becomes 

At high-temperature the spin-drag relaxation rate be- 
comes C-independent and identical to the one calculated 
in Ref. [12 for a strictly unpolarized system. 



input for the application of time-dependent current-spin 
density functional theoryi^ to the study of the propaga- 
tion of spin pulses in the nonlinear regime^^. This in turn 
will be useful in future comparisons between theory and 
experiments on spin-pulse propagation in ID cold Fermi 
gases. 

Note added in proof — Going beyond second-order per- 
turbation theory, one finds to leading order an additional 
contribution to 1 /Tgd that arises from intra-spccics inter- 
actions. Such a contribution has a power-law tempera- 
ture dependence^Si^, and could thus become dominant 
over the exponential contribution calculated in this work 
at the lowest temperatures. In Ref. [13 it has been shown 
that this contribution goes as , and is proportional to 
^o(^o — i'2fcF)^- This implies that if the interaction range 
is less than the typical interparticle distance the fourth- 
order contribution is expected to be suppressed. 



IV. CONCLUSIONS 

In summary, we have carefully studied within second- 
order perturbation theory the spin-drag relaxation time 
of a spin-polarized ID Fermi gas. We obtain accurate nu- 
merical results for 1 / Tgd as a function of spin polarization 
and temperature (see Fig. [T]), and also an accurate an- 
alytical formula for the spin-polarization dependence of 
1/Tsd in the low temperature regime, which is completely 
different from what was obtained earlier from an approx- 
imate argumcnt*2£. These results provide the necessary 
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